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Abstract 

Using the long wave perturbation scheme (gradient expansion), the effect 
of inhomogeneity on the inflationary phase is investigated. We solved the 
perturbation equation of which source term comes from inhomogeneity of a 
scalar field and a seed metric. The result indicates that sub-horizon scale 
inhomogeneity strongly affects the onset of inflation. 



C^; I. INTRODUCTION 

52 I Dynamics of inhomogeneity plays an important role in many phenomena in general rela- 

^ tivity. Generality of inflation is one of such a problem. To explain the present our universe, 

"X' 

Vh ' we must have an inflationary phase in past. But as we do not know the initial condition 

of universe, we hope that universe can enter inflationary phase from wide range of initial 
condition. For homogeneous closed universe, the universe cannot enter inflationary phase if 
the initial curvature is too large [1,2]. For inhomogeneous universe, the effect of strong (non- 
linear) initial inhomogeneity is studied for space-time with symmetry by numerical method 
[7]. These result indicate that sub-horizon scale inhomogeneity strongly prevent the on- 
set of inflationary phase. But imposing symmetry to space-time may lose general feature 
of dynamics of inhomogeneity. We want to know about more general situation without 
symmetry using appropriate approximation that include non-linear effect. Gradient expan- 
sion(GE) [3-6] is one of such a method. This approximation scheme expands the Einstein 



equation by the number of the spatial gradient. As the background solution, we have an 
inhomogeneous space-time that does not include the effect of a spatial curvature. We can 
include the curvature by calculating next order. This method describes a long-wave non- 
linear perturbation and it is suitable to applying to the inflationary universe. Furthermore, 
we can treat non-linear effect without imposing any symmetry for a space-time. 

In this paper, we treat Einstein gravity with a minimally coupled scalar field by GE. Our 
purpose is to observe the effect of large scale inhomogeneity on inflationary phase. Especially, 
we are interested in how the onset of inflation is affected by the initial inhomogeneity. 

This paper is organized as follows. In Sec. II, we shortly review GE approximation using 
Hamilton- Jacobi method. We derive the second order equation of motion that include the 
scalar field inhomogeneity. In Sec. Ill, we obtain the solution for the inflationary phase and 
discuss the effect of inhomogeneity on inflation. Sec. IV is devoted to summary. 

II. EQUATION OF MOTION BY H-J METHOD 

We follow the method of Salopek [5] and derive the second order equation of motion of 
GE by Hamilton- Jacobi(H- J) method. H-J equation is 

hijhki 2 2k 2 dcj) 

5S 5S 5S . , 

Hi = -2{'jik- — )j + - — 7«fc,j + Y-0,i = 0, (1) 

where k = 8ttG. These two equations are Hamiltonian constraint and momentum constraint, 
respectively. We expand the generating functional S by the number of spatial derivative 

The form of generating functional S^'^''^^ is assumed to be 

5(°) = -- / d?x-i^l^H{cP), 5(2) = j d'x-i^'^{J{(P)R + K{(j))(f).i(f)''), (2) 



where ; represents spatial covariant derivative associated with 7jj. This form of generating 
functional satisfies momentum constraint strongly provided that parameter fields (integration 
constant of H-J equation) are spatially constant. 
The zeroth order Hamiltonian becomes 
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and the requirement that it vanishes gives 
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The second order Hamiltonian becomes 
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The requirement that this vanishes gives the following equations for J, K: 
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HJ 



2k k 



-H J =0, 



-H'K - Hj' = 0, 

K 

HK + - + -H'k' - AHJ" 

2 K 



(5) 



0, 



where = ^. We have two equations for K^ but it can be shown that the second equation 
automatically satisfies the third one. So they are consistent. 

Now the evolution equation that is accurate to the second order spatial gradient becomes 
the following form. From now on, we use the synchronous coordinate system N = l^Ni = {): 

,6S 
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The first order solution is obtained by integrating the following equation: 



-V(0), 

K 



7ij = 2H-iij, 



(8) 



and the solution becomes 

where hij{x) is an arbitrary spatial function(seed metric) and A^^'^ = J dtH. to is a inte- 
gration constant that depends only on spatial coordinate, and a scale factor A depends on 
spatial coordinate through 0. We can rewrite above equations (6), (7) using to and covariant 
derivative associated with a seed metric hij. The spatial curvature becomes 



i?,,(7(°)) = R^J{h) - h,, [{H + //')to|„tol™ - HAt 
+ {-H + H )to|jto|j + HtQ\ij, 



(9) 
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i?(7(°)) = -{R{h) + {-AH - 2ij2)t t^l™ + 4iJAto), 



where tou, means covariant derivative associate with hij. The derivative of 0, J is 
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To derive these expression, we have used eq.(5). Using these relations, the form of evolution 
equation becomes 



0=--iJ' + ^i?(/.)-i:^to|„tol™, 
'jij = 2H-fij + AKJ{-R{h)hij - Rij{h)) 

+ 2H [—tQ^itQ^j + -hijto\mto ) — 2to|ij. 



(12) 
(13) 



We consider the perturbation due to the second order spatial gradient. By linearizing 
eq.(12) and (13), we get 



S^ = --H"6(P + E{(Po,h), (14) 

K, 

5Kj=2H'hi,5(t) + Fi,{(t)o,h), (15) 

where 5hij = 6'~fij/A and E, Fij are source terms from the second order gradient: 

7/ 1 TTl 

E = —rR{h) ^^oim^o , 

A K, A ' 

7 1 O" 1 O 

Fjj = AK,— {-R{h)hij — RijQi)) + 2— (— to|iio|j + ■7:hijtQ\m^Q ) — ^^o|jj- 

These terms represent non-hnear inhomogeneity due to the seed metric hij and the zeroth 
order matter field (f){t — to (a;)). The solution of (14) becomes 

* E\ ■ 
dt^ 00, (16) 

ti 00/ 

where we have chosen the integration constant such as 6(f){ti) = 0. Using this solution, the 
metric perturbation becomes 



Shij = f {2H'5ct>hj + F,j) 

Jti 

= 2H dt^hij + dti -2H^hij + Fij . (17) 

Jti 00 "'** V 00 / 

III. SOLUTION IN THE INFLATIONARY PHASE 

We derive the solution in the inflationary phase a > 0. We flrst evaluate the function 
J: 

J = ( / dta + const ) 

2Ka \J J 

= —{H~^ - H~\H-'y + h-Hh-\h-^\\ 

2k 

- \ {a^Hr^ - a,Hr\Hr^y + ...) + ^}, (18) 

and we choose the integral constant so that the term proportional to 1/a vanishes. 

J = i-ij-i {i - {H-'y + {H-\H-'yy -•■■}. (19) 



During inflation, we can neglect all but the first term in this expression because the condition 
|i/| < H^ is satisfied. To the leading order of this expansion, we can evaluate the following 
integrals: 
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and the perturbation of the scalar field becomes 
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(20) 



Therefore during infiation, 
Hubble time scale: 



approaches a constant value and 50 goes to zero in a 
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The metric perturbation during infiation is expressed as 
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(21) 
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The metric that is not proportional to hij approaches to constant values in a Hubble time 
scale. GE is correct if |0o| > \S(f)\ and \hij\ > \Shij\ are satisfied. These conditions are 
rewritten to be 



R{h) toi^tpl'" 



2am^ 



<2, 



|i?(7(°))| < \RA = 12H\ 



(23) 



They are well satisfied during infiation provided that they are satisfied at onset of infiation 
because the quantities in the left hand side of the condition decrease in time during infiation. 
Therefore the expression of perturbation (21), (22) are correct in whole time of infiation. 

A. Effect of the spatial curvature on inflation 

We investigate the effect of the curvature inhomogeneity. For this purpose, let us consider 
a local Hubble radius and a local scale factor that are defined by 
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a = exp( / dtH). 



(24) 



Using eq.(12), these quantity become 

H = h(i 
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The co-moving Hubble radius becomes 
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This quantity decreases during inflation and grows in non-inflationary phase. We choose 
ti,tf are onset and end time of inflation, respectively. At t = tj, {aH)~^ has maximum and 
at t = tf minimum. The maximum value at t = ti for R{x ) < becomes smaller than one 
for R('j^'^^) = and larger for R('j'^'^^) > 0. This implies that the negative -R(7*-°-') enhances 
and the positive -R(7*^°^) suppresses inflatinary phase. Especially if i?(7'^°)) becomes as large 
as -R* = 12-?/^^j, the co-moving Hubble radius goes to infinity at t = tj and inflation may 
not occur. Of course we cannot apply GE to such a regime because the condition (23) is 
violated. But we extrapolate the behavior for small \R\ and predict the behavior for large 
\R\. As the curvature radius is deflned by Lr - 



'6/R, this critical value corresponds to 
the scale Lr = \/2H^^. If the scale of curvature inhomogeneity is larger than this scale, 
inflation may occur. But the scale is smaller than this value, we may not have inflation. 
This result is consistent with the result of the calculation of numerical relativity [7]. 

To see the effect of curvature from a different view point, we consider local e-folding time 
of inflation: 



N{U 



ti 



tf - 1 

dtH ^ N ^ 

ti 24 
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H^ 



tf 



N - 
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24' H"^ 



(27) 



where N is the e-folding for R{x ) = 0. Positive R makes e-folding shorter and negative 
R makes it longer than R = case. This also indicates that positive spatial curvature 
suppresses inflationary phase. 
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We check the strong energy condition: 



V- 6 ^ . / H ,^ R(h) 



U\mtr) . (28) 



12(a2if2) ■ 602°!"^° J 

At the onset of inflation, {H / H'^)ti = -1, so (p + Y.Pi)ti = H^ / K{R{h) / {20^ H'^)ti + 
^oim^o'^/la^)**)- Positive R{h) or to|m^o'™ term make (p + Y.Pi) > and the strong en- 
ergy condition is satisfied. This means that inflation does not occur. Negative R{h) can 
make (p + '}2Pi)ti < if its value is larger than the to|m.^o''" term, and the strong energy 
condition is violated. 

B. Numerical calculation 

To check the result of the previous subsection, we solve eq.(13),(14) numerically. We 
prepare the initial condition at pre-inflationary phase. For simplicity, we consider only 
to = case. The model is V"(0) = l/2m?(jp. Parameter and initial value is m? = 0.5, 0(0) = 
7, 0(0) = —5, J(0) = in K = 1 unit. The scalar curvature R{h) = —30, 0, 30. This value 
is chosen to satisfy the condition (23). Fig.l is a phase space diagram (0,0). The solid 
line is R{h) = 0, dotted line is R{h) < and dashed line is R{h) > 0. The zeroth order 
solution (i?(/i) = 0) has inflationary phase from ^ 6 to ~ 1. For positive R{h), the 
trajectory enters inflationary phase with smaller value(smaller H). For negative R{h), 
the trajectory enters inflation with larger value(larger H). After entering inflation, all 
trajectories approaches attractor solution. Fig. 2 is time evolution of the perturbation of 
scalar fleld (50. In Hubble time scale, its value approaches a constant value. Fig. 3 shows the 
co-moving Hubble radius. R{h) = curve has a maximum at t ~ 0.06 and after this time, 
the system enters inflation. The positive R{h) makes the co-moving Hubble radius larger 
and the negative R{h) makes smaller. These behavior is consistent with the analysis of the 
previous subsection. 



IV. SUMMARY 

As we have shown, GE becomes good approximation in the inflationary phase because the 
perturbation goes to a constant value. If the initial inhomogeneity is not so large as to break 
the approximation, we can use GE as a model to construct an inhomogeneous inflationary 
space-time. As is well known, a global structure of inflation is very complicated because 
of continuous creation of new inflationary domain via quantum fluctuation. We can grasp 
these feature by using a stochastic approach [8-10], but this method does not incorporate 
the information of geometry completely. GE can describe the space-time of inflation whose 
local Hubble radius is different from place to place. The zeroth order solution has the form 
(f) = (f)Q{t~tQ{x)) and choosing the time delay function to{x) appropriately, we can construct 
such a space-time as the solution of GE. We will discuss on this topics in a separated 
publication [11]. 



Note: After this work has completed, we noticed that the paper [12]. They treat this 
subject in the same context as ours. But they perform only numerical analysis and does not 
derive analytic expression of perturbation. 
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FIGURES 
FIG. 1. A phase space diagram (<^, </>). for R{h) = O(solid line), R{h) > O(dotted line) and 
R{h) < 0(dashed line). R{h) = solution enters inflation at </> ~ 6. For positive R{h), solution 
enters inflation with smaller value and for negative R{h), solution enters inflation with larger (p 
value. 

FIG. 2. Time evolution of perturbation (5</>. Perturbation approaches constant value in Hubble 
time scale. 

FIG. 3. Time evolution of co-moving Hubble radius. System enters inflationary phase at 
t ~ 0.06 for R{h) = case. For positive R{h), co-moving Hubble radius at the onset of infla- 
tion becomes larger. For negative R{h), co- moving Hubble radius at the onset of inflation becomes 
smaller. 



11 



vphi 



-0.45 



-0.5 



-0.55 



-0.6 



-0.65 




phi 



Fig.l 



dphi 





/ 


-^ ■ 


^_ ^^^ ^^^ 








. 6 








"— ^ '-—m. 


0.4 


/ 














/ 












0.2 


/ 


















2 


4 


6 


8 10 


0.2 














0.4 


- \ 
\ 














\ 












n r 


\ 
\ 
\ 










___--" 


. 6 


1* ^ 











Fig. 2 



1/ (aH) 



0.35 



0.3 



0.25 



0.2 



0.15 




Fig. 3 



